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Atomic hydrogen provides a unique test case for computational electronic structure methods, since 
its electronic excitation energies are known analytically. With only one electron, hydrogen contains 
no electronic correlation and is therefore particularly susceptible to spurious self-interaction errors 
introduced by certain computational methods. In this paper we focus on many-body perturbation- 
theory (MBPT) in Hedin's GW approximation. While the Hartree-Fock and the exact MBPT 
self-energy are free of self-interaction, the correlation part of the GW self-energy does not have this 
property. Here we use atomic hydrogen as a benchmark system for GW and show that the self- 
interaction part of the GW self-energy, while non-zero, is small. The effect of calculating the GW 
self-energy from exact wavefunctions and eigenvalues, as distinct from those from the local-density 
approximation, is also illuminating. 

PACS numbers: 31.25.Jf, 31.15.Lc, 31.15.Ar 



I. INTRODUCTION 

Ab initio many-body quantum mechanical calculations 
are crucially important to our understanding of the be- 
havior of atomic, molecular and condensed matter sys- 
tems. It is well known that prediction of the behavior 
of these systems requires the description of electronic 
correlation. Whilst density-functional theory (DFT) in 
the local-density approximation (LDA) does this with 
startling success in many cases, it does so at the expense 
of a non-physical electron self-interaction. For delocal- 
ized electron systems this self-interaction becomes negli- 
gible, but in atomic or strongly localized electronic sys- 
tems it plays an important role. If one is interested in 
the calculation of quasiparticle excitation spectra, many- 
body perturbation-theory (MBPT) is formally a correct 
way to proceed. For solids, MBPT in Hedin's GW ap- 
proximation [l[ has become the method of choice, but it is 
also increasingly being applied to molecular systems and 
clusters. The GW self-energy can be decomposed into 
correlation and exchange parts, where the latter is the 
same as the Fock operator encountered in Hartree-Fock 
theory and thus self-interaction free. While the exact 
self-energy must also be free of self-interaction, the cor- 
relation part of the GW self-energy does not have this 
property. To investigate the influence of self-interaction 
in the GW approach the hydrogen atom provides an ideal 
case, because the exact solution is known analytically. 

Hydrogen in its solid phase has previously been stud- 
ied within the GW approximation by Li et al. Q, who 
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analyzed the transition between the high-pressure solid 
phase and the low density, atomic-like limit. For indi- 
vidual atoms, GW electron removal and addition ener- 
gies (we use the term "quasiparticle" energies by analogy 
with the solid-state situation) have been investigated by 
Shirley and Martin H, Dahlen et al. 0,11], Stan et al. Q 
and Delaney et al. [7[, although hydrogen was not con- 
sidered. These studies have shown that GW, in general, 
gives quasiparticle properties which are much improved 
over DFT and Hartree-Fock methods, even for atoms. 

In this paper we use the hydrogen atom as a benchmark 
system to quantify the self-interaction error in the GW 
approach. Since the self-energy diagrams beyond GW, 
known as the vertex correction, must by definition cor- 
rect this self-interaction error, our findings are relevant 
for research into vertex functions for the many-electron 
problem. 

Attention has recently focused on the prospects for im- 
proving the usual non-self-consistent GW calculations by 
choosing an initial Green's function, Go, that is physi- 
cally more reasonable than the LDA (e.g. @,i,[l). We 
explore this here by determining the sensitivity of the 
self-interaction error to the use of the exact hydrogenic 
orbitals and energies in place of those from the local- 
density approximation (LDA). We also assess the error 
introduced into GW calculations by employing first-order 
perturbation theory in solving the quasiparticle equation 
(as opposed to the full numerical solution), and we ana- 
lyze the quasiparticle wavefunctions that emerge from a 
full solution. 



II. HARTREE-FOCK VS. DFT-LDA 

In many-body perturbation theory the quasiparticle 
excitation energies and wavefunctions ipi a are the so- 
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lutions of the quasiparticle equation 

a' 

(i) 

where, in Hartrcc atomic units, i?o( r ) = — 5 + v ext(?) 
and v ext (r) is the external potential. It is customary 
to divide the mass operator M into the local Hartree 
potential (vh) and the non-local self-energy (E) 

AO(r, r'; e) = v H (r)S(r - r')S aa , + E CTCT , (r, r'; e). (2) 

Omitting correlation contributions from E yields the 
exact-exchange or Hartrce-Fock case, where the self- 
energy takes the form 
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and the Hartree potential is given by 
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Since the sum runs over all occupied states the Hartree 
potential contains an artificial interaction of electrons 
with themselves. For the hydrogen atom this so-called 
self-interaction of the electron is the only content of the 
Hartree potential and may be calculated analytically us- 
ing the exact ground-state wavefunction 
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potential v xc , which is formally given as the functional 
derivative of the exchange-correlation energy E xc : 
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In analogy to the quasiparticle equation |T]) the Kohn- 
Sham eigenvalues e^. s are often interpreted as excitation 
energies, although this is not formally justified. 

One of the most common approximations for E xc is 
the local-density approximation (LDA) [l(| , in which the 
many-body exchange and correlation contributions to the 
total energy (E xc ) are included by comparison with the 
homogeneous electron gas (HEG): 
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Here \12\ w e follow the parameterization of Perdew and 
Zunger 13| for the exchange-correlation energy density 
e™ G [n(r)] of the homogeneous electron gas based on the 
data of Ceperley and Alder [3] . 

The LDA has been remarkably successful at accounting 
for correlation even in systems that are highly inhomoge- 
neous. However, it is well known that the introduction of 
correlation in the LDA comes at the expense of the exact 
treatment of the self-interaction. Because the exchange 
functional is taken from the homogeneous electron gas 
it no-longer cancels the spurious self-interaction present 
in the Hartree term. In most systems this is a minor 
effect and is more than compensated by the improved 
treatment of the electron correlation. The LDA can be 
improved by explicitly removing the self- interaction [l3| , 
however this becomes increasingly difficult as the sys- 
tem's complexity is increased. 
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The self-interaction is a positive, decreasing function and 
hence tends to delocalize the wavefunction (i.e. incor- 
rectly pushes its weight away from the nucleus). 

In Hartree-Fock the self-interaction terms introduced 
in the Hartree potential are exactly canceled by the Fock 
operator E^. This makes Hartree-Fock exact for one- 
electron systems such as the hydrogen atom. However, 
the lack of correlation renders Hartree-Fock unsuitable 
for many polyatomic systems of interest. 

Of more practical use is Kohn-Sham density- 
functional theory, which by virtue of the Hohenberg- 
Kohn theorem [1JJ establishes an exact and universal 
functional relationship between the ground-state density 
n(r) and the total energy E of a system. Mapping the 
system of interacting electrons onto a fictitious system 
of non-interacting electrons, that reproduces the exact 
density, yields the Kohn-Sham equations: 

[H (t) + v H (v) + v xc (r)} <Mr) - <& £ WW- (?) 

All electron-electron interactions beyond the Hartree 
mean field are encompassed by the exchange-correlation 



III. THE GW APPROXIMATION 

In Hcdin's GW approximation [l[ the self-energy in 
Eq. [2] is given by 
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where S is an infinitesimal positive time. At the level of 
GW, spin flips are not accounted for and the input 
Green's function is diagonal in its spin representation 
G(j(j — Gap/ 5(j(j' 

In common with usual GW calculations, a Kohn-Sham 
Green's function G° is used for G, given by 
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and makes the non-interacting polarizability 
X° CT (r, r', e) = — / de'G° CT (r, r'; e' - e)G° ff (r', r; e') 
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and the dielectric function spin-dependent (though spin- 
diagonal) . The inverse dielectric function in the random- 
phase approximation 



S(r-r')- f dr%(r,r")5>°a(rV',e) 

(13) 



and thus the screened Coulomb interaction 



W (r,r',e)= f dr"e-\r, r", e)v(r", r'). (14) 

then emerge as spin- independent quantities, giving rise 
to the simple spin dependence in the GW self-energy 
(Eq.HnD. 

For numerical convenience and physical insight we sep- 
arate the GW self-energy (jlOp according to 
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S CTCT = -i[G„v + G aa (W - «)] = + K 



The first term (E* a ) corresponds to the Fock operator 
in Eq.[3]and will exactly cancel the self-interaction intro- 
duced by the Hartree potential. It is therefore immedi- 
ately clear that any deviation from the exact result for 
hydrogen can only come from the correlation part of the 
self-energy (££ ff ). 

The incorrect self-interaction affects the electron re- 
moval energies (here, the ionization potential). For elec- 
tron addition energies such as the electron affinity, the 
entire Hartree potential has a physically reasonable inter- 
pretation, since it acts on the wavefunction of the origi- 
nally unoccupied state which has not contributed to the 
electron density. 



IV. COMPUTATIONAL APPROACH 

We solve the quasiparticle equation JTJ) with the GqWq 
self-energy (fT0|) for the quasiparticle energies and wave- 
functions by fully diagonalizing the quasiparticle Hamil- 
tonian in the basis of the single particle orbitals of the 
non-interacting system. Since the ground state of the 
hydrogen atom (Eq. [5]) is spherically symmetric, it is 
sufficient to describe all non-local operators in the GW 
formalism by two radial and one spin coordinates, r,r',a 
and one angular coordinate, 9, that denotes the angle 
between the vectors r and r'. The self-energy (Eq. [T0|) 
then assumes the much simpler form 
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where P/(cos6*) is a Legendre polynomial of order I. 

The Legendre expansion coefficients of the self-energy 
are calculated directly, thereby surpassing the need for a 
numerical treatment of the angular dependence. We use 
a real-space and imaginary time representation to 
calculate the self-energy from the non-interacting Green's 
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FIG. 1: Radial probability distributions of the hydrogen Is 
state: the quasiparticle wavefunction deviates only slightly 
from the exact wavefunction, when the latter is used as a 
starting point (exact+GIT). The LDA wavefunction, on the 
other hand, is more delocalized as a result of the inherent self- 
interaction. Adding quasiparticle corrections (LDA+G1T) 
brings the resulting quasiparticle wavefunction slightly closer 
to the exact one again. The inset shows the difference to the 
exact wavefunction. 



function Go- The expression for the self-energy on the 
real frequency axis is obtained by analytic continuation 
fl6( |. The current implementation has been successfully 
applied to jellium clusters [l~7| and light atoms 

Our code allows us to solve the quasiparticle equation 
(Q]) for the GW self-energy with no further approxima- 
tion. However, in order to separate the contribution that 
arises from the correlation part of the self-energy from 
that of the exchange part and the Hartree and exchange- 
correlation potential we also solve the quasiparticle equa- 
tion with the frequently made approximation that the 
quasiparticle wavefunctions are given by the Kohn-Sham 
wavefunctions. The resulting equation for the quasipar- 
ticle energies is 
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where the brackets () denote matrix elements with re- 
spect to the Kohn-Sham wavefunction <^j CT . 

In order to explore the role of the starting points for a 
GW calculation, two possible Kohn-Sham input Green's 
functions are chosen. First, the familiar LDA, and, sec- 
ond, the exact Kohn-Sham solution for the hydrogen 
atom which has the exact wavefunction of the hydrogen 
Is state ([5]) and v xc (r) = —vh(v). (This exact Kohn- 
Sham Green's function, incidentally differs from the ex- 
act Green's function of the hydrogen atom, because the 
exact Kohn-Sham unoccupied eigenvalues do not signify 
electron affinities. The exact Green's function cannot 
be constructed from any orthonormal set of one-particle 
wavefunctions. ) 
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Exact 


HF 


LDA 


LDA+GW 


Exact+GW 


-13.61 


-13.61 


-6.36 


-12.66 


-13.40 



TABLE I: Quasiparticle energies (eV) for the Is state of hy- 
drogen (the ionization potential) obtained by diagonalizing 
the quasiparticle Hamiltonian {T}. Two GW calculations are 
shown, starting from the LDA and from exact Kohn-Sham, re- 
spectively. For comparison, the Hartree-Fock (HF) and LDA 
eigenvalues are also shown. 



Kohn-Sham Go 


GW 




(Vxc) 


<£c> 


LDA 


-12.93 


-15.38 


-8.25 


0.56 


Exact 


-13.35 


-17.00 


-17.00 


0.25 



TABLE II: Quasiparticle energies (eV) for the Is state of hy- 
drogen obtained by solving Eq. (|17[) . The contributions from 
the exchange (E x ) and correlation (E c ) part of the self-energy 
are compared to that of the exchange-correlation potential 
(vxc) for the LDA and the exact case (v xc = —vh) as a start- 
ing point. Exact value for ei s is —13.61 eV. 

V. RESULTS AND DISCUSSION 

The calculated ionization potentials (from a full solu- 
tion of the quasiparticle equation) are shown in Table 
HJ The self-interaction errors in the two GW quasipar- 
ticle energies are seen to be fairly small: 0.95 eV when 
the approximate LDA Kohn-Sham starting point is used, 
and the much smaller 0.21 eV when the exact Kohn- 
Sham starting point is used. Clearly the LDA is such a 
physically poor representation of the correct physics in 
this extreme system (owing to the large self-interaction 
present in the LDA calculation itself, as reflected in the 
large error in the LDA Kohn-Sham eigenvalue) that it 
forms a very unsuitable starting point for GW. How- 
ever, a physically reasonable starting point reduces the 
GW self-interaction error to a small size. 

Since (S c ) gives a non-vanishing contribution to the 
hydrogen Is state, even if the analytic solution is used as 
a starting point, the quasiparticle wavefunction will differ 
from the exact one. Figure [1] shows that the GW corre- 
lation gives rise to a slight derealization of the quasipar- 
ticle wavefunction in this case. This relaxation, however, 
now makes the quasiparticle wavefunction an eigenfunc- 
tion of the quasiparticle Hamiltonian and reduces the 
deviation from the exact energy of the Is state to 0.2 eV, 
as shown in Tabic HI In the LDA the self- interaction error 
is much more pronounced and the wavefunction becomes 



significantly more dclocalized. The GW self-energy cor- 
rects this to a small extent (as reflected in the quasipar- 
ticle wavefunction), but the remaining discrepancy reit- 
erates the unsuitability of the LDA as a starting point 
for GW in this self-interaction-dominated atom. 

For an analysis of the contributions to the self-energy 
we turn to the perturbative solution of the quasiparti- 
cle equation using Eq. [T71 shown in Table fll] When 
the exact Kohn-Sham wavefunction and eigenvalues are 
used, as in the Hartree-Fock case the exchange part of the 
self-energy is seen to cancel the self-interaction contribu- 
tion from the Hartree potential exactly. The correlation 
part, on the other hand, is not zero, but amounts to a 
self-polarization of 0.25 eV. When the LDA is used as 
starting point the influence of the LDA wavefunction on 
the exchange operator becomes apparent and it reduces 
from —17.00 eV in the exact case to —15.38 eV. This cor- 
rects the highly overestimated LDA eigenvalue for the Is 
state of -6.36 eV (see Tabled to -13.49 eV. However, in 
this case the contribution from the correlation part of the 
GW self-energy is even larger than when starting from 
the exact case and increases the quasiparticle energy to 
-12.93 eV. 



VI. CONCLUSION 

We have performed spin-resolved benchmark calcula- 
tions for the GW formalism using the analytically known 
solutions of the hydrogen atom as a reference, making the 
self-interaction error introduced by the correlation part 
of the GW self-energy directly assessable. When the ex- 
act Kohn-Sham Green's function is used as the input 
to GW, the self-interaction error is small (0.21 cV, one- 
thirtieth the size of that in the LDA), but not negligible. 
If the LDA Kohn-Sham Green's function is used, as done 
in many GW calculations for more complex systems, a 
larger self-interaction error remains, inherited from the 
LDA starting point. 
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